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摘 要 : 循环 荷载 作用 下 深厚 饱和 软 番 士 的 固 结 与 长 期 沉降 变形 分 析 对 公路 、 地 铁 等 基础 设施 的 沉降 分 析 有 很 重要 的 
意义 。 传 统 的 固 结 分 析 方 法 无 法 合理 地 考虑 到 长 期 循环 荷载 作用 下 软 兢 土 的 蠕 变 沉降 ， 而 基于 弹 粘 塑 性 的 本 构 模型 分 
析 又 不 能 考虑 长 期 循环 荷载 作用 下 土 体 固 结 过 程 中 渗透 系数 ky 和 体积 压缩 系数 m 的 非 线 性 变化 。 本文 将 综合 考虑 长 期 
循环 荷载 作用 下 饱和 软 黏 士 的 非 线 性 固 结 及 其 固 结 过 程 中 产生 的 蠕 变 。 本 文 构建 了 一 种 具有 普遍 适用 性 的 一 维 弹 黏 塑 
性 非 线 性 固 结 模型 ， 可 以 较 好 的 计算 长 期 循环 荷载 作用 下 考虑 蠕 变 效应 的 饱和 软 秋 士 的 长 期 沉降 变形 。 本 文 与 土工 试 
验 相 结合 , 给 出 了 便于 应 用 的 模型 参数 确定 方法 ; 并 提出 了 一 种 多 级 加 载 过 程 中 的 初始 条 件 确定 方法 , 该 方法 可 以 反映 
出 应 变 和 超 静 孔 压 在 深度 上 的 不 均匀 性 。 采 用 杭州 软 秋 士 进 行 了 多 组 不 同 振幅 不 同 频率 不 同 动态 偏 移 量 的 循环 荷载 固 
结 试验 ， 通 过 计算 结果 与 试验 结果 的 对 比分 析 ， 验 证 了 模型 的 合理 性 。 
关键 词 : 循环 荷载 ， 蜂 变 ， 非 线性 固 结 ， 固 结 试验 ;沉降 计算 


One dimensional elastic visco-plastic nonlinear consolidation model of soft clay 
under cyclic loading 


JIA Dongge, Zheng Yangguang 
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Abstract: Long-term consolidation and settlement analysis of soft clay under cyclic loading is very important for settlement 
prediction of infrastructures such as highways and railways. The traditional consolidation analysis method does not take into 
account of creep settlements. However, the constitutive model for creep analysis will not suit for the consolidation analysis under 
cyclic loading with nonlinear behavior of permeability coefficient kv and volume compressibility coefficient mv. This study 
establishes a one-dimensional elastic visco-plastic nonlinear consolidation model with universal applicability for cyclic loading 
conditions. The proposed model could yield an accuracy result which considers the long-term creep and nonlinear behavior of 
soft soil. The determination of model parameters is also analyzed. In addition, a method of determining the initial conditions in 
the multistage loading process is proposed by polynomial fitting. A series of one-dimensional cyclic loading tests are conducted 
in this study and verified with the proposed model. Comparisons reveal that the proposed method can reflect the inhhomogeneity 
of the strain and the excess pore pressure in the depth. The model could also capture the consolidation behavior of soft clay under 
different loading amplitudes, frequencies and dynamic offsets. The rationality and effectively of the proposed model is discussed 
and verified with the experimental results. 


Key words: cyclic loading; creep; nonlinear consolidation; consolidation experiment; settlement calculation. 


0 引 站 由 忽略 里 变 特性 的 固 结 分 析 得 出 的 土 体 沉降 会 明显 小 
于 实际 情况 ， 影 响 工程 的 安全 性 。 

循环 荷载 是 工程 中 常见 的 一 种 荷载 形式 ， 波 浪 、 上 个 世纪 60 年 代 以 来 ， 不 少 学 者 开始 对 土 的 蜂 

车 辆 、 粮 仓 、 地 下 水 等 周期 性 变化 都 会 引起 土 层 承受 ” 变 特 性 开展 研究 。Walker、Mesri 等 (8-10] 通 过 对 室内 试 


周期 性 荷载 。 已 有 不 少 学 者 忽略 土 的 非 线性 ， 对 此 情 


由 


况 下 的 固 结 特性 进行 了 研究 t-g。 本 世纪 初 , 谢 康 和 加 


验 数据 的 分 析 得 出 了 粘土 的 蜂 变 率 与 压缩 指数 近似 成 
线性 关系 的 结论 。 息 建华 建立 了 蠕 变 参数 与 时 间 的 关 


在 土 体 渗透 性 的 降低 与 压缩 性 的 减 小 成 正比 的 假设 下 ，“ 系 01， 克服 了 肾 变 参数 恒定 时 里 变 将 随时 间 的 增 大 而 
给 出 了 循环 荷载 作用 下 地 基 一 维 非 线性 固 结 的 解析 解 ，“ 无 限 增 大 的 不 足 。Singh 等 建立 的 蠕 变 模型 通过 调整 
但 以 上 研究 均 没 有 考虑 土 的 蠕 变 特性 。 对 于 软 黏 十 ， 参数 m 与 1 的 相对 大 小 可 以 分 别 表现 出 里 变 率 与 时 间 


对 数 的 负 相 关 、 无 关 和 正 相 关 关 系 02]。 用 于 描述 士 的 
需 变 特性 的 流 变 模型 也 得 到 了 一 定 的 发 展 ， 主 要 有 马 
克 斯 威 尔 模 型 ， 开 尔 文 模型 和 宾 厄 姆 模型 3-5。 早 在 
1967 年 ，Bjerrum 提出 了 一 种 等 时 间 线 模型 来 解释 土 
的 延迟 压缩 0 。 之 后 ，Garlanger 基于 等 时 间 线 模型 给 
出 了 特定 应 力 应 变 情 形 下 土 的 应 变 计算 方法 01。 括 建 
华 和 Graham 对 等 效 时 间 进 行 了 物理 解释 与 数学 定义 ， 
导出 了 具有 普遍 适用 性 的 一 维 弹 黏 塑性 (EVP) 本 
模型 0829。 为 了 考虑 士 体 变形 中 的 固 结 特性 ， 括 建 
华 等 进一步 运用 有 限 差 分 和 有 限 单元 法 将 一 维 EVP 
模型 与 一 维 固 结 进行 了 耦合 分 析 P2。 
维 弹 黏 塑性 本 构 模 型 合理 的 考虑 了 士 的 流 变 性 ， 
且 具 有 一 定 的 普遍 适用 性 。 而 近 些 年 来 ， 对 循环 荷载 
作用 下 土 体 的 固 结 理论 多 忽略 土 体 的 蠕 变 。 本 文通 过 
将 非 线 性 固 结 理论 与 EVP 本 构 模 型 进行 耦合 , 得 出 了 
考虑 土 体 流 变 特性 的 一 维 非 线 性 固 结 模 型 。 之 后 对 杭 
州 饱和 软 黏 士 进行 了 多 组 一 维 循环 荷载 固 结 试验 ， 将 
试验 结果 与 用 Mathematica 计算 的 模型 结果 进行 对 比 
分 析 ， 验 证 了 循环 荷载 作用 下 一 维 弹 莫 塑性 非 线 性 固 
结 模型 的 适用 性 。 


时 豌 


1 模型 的 建立 


1.1 一 维 非 线 性 固 结 理论 

图 1 为 本 文 拟 求解 的 计算 模型 简 图 。 厚度 为 瓦 的 
饱和 软 厅 土 受 均 布 荷 载 作 用 ，p(?) 为 任意 形式 的 循环 
荷载 。 考 虑 到 土 在 固 结 过 程 中 的 渗透 性 是 变化 的 ， 采 
用 式 (1) 计算 土 体 在 不 同 应 力 应 变 状 态 下 的 渗透 系数 


D2]。 


1 和 荷载 及 边 珊 情况 


Fig. 1. Loading and boundary condition of soil layer 


二 I (1) 
Oo; 


式 中 : 06、kw 为 e 一 1go, 和 elgk, 曲线 上 任意 


“0o” 点 处 的 有 效应 力 和 渗透 系数 ，cc 为 压缩 指数 ， ck 


为 渗透 指数 。 
饱和 土 体 的 一 维 固 结 连续 方程 为 
0 | k, Ou Oe, 
v = 2 
引 yy | ot ® 


将 式 (1) 代入 式 (2) 得 : 


求 导 展开 为 
ko | cf 1 YouV| ae 
| [oe ee 
式 (4) 即 为 考虑 土 固 结 过 程 中 的 渗透 系数 和 压缩 
系数 非 线 性 变化 而 得 到 的 控制 方程 ， 其 中 o'=p-u。 


1.2 ”一 维 弹 黏 塑性 本 构 模 型 
笑 建 华 和 Graham 提出 的 一 维 弹 黏 塑性 本 构 模 型 

是 一 种 非 线 性 流 变 模型 0&829。 如 图 2 所 示 ， 在 土 的 应 
力 应 变 关系 中 采用 了 “有 瞬时 压缩 时 间 线 ” “参考 时 间 
线 ” “等 效 时 间 线 ”等 概念 。 土 体 沿 着 瞬时 压缩 时 间 
线 产 生 的 应 变 是 完全 弹性 的 ,不 计 入 时 间 因 素 的 影 
不 产生 蠕 变 。 等 效 时 间 线 是 一 组 平行 线 ， 每 一 条 线 代 
表 相 同 的 蠕 变 率 。 蠕 变 率 的 大 小 可 以 用 等 效 时 间 万 表 
示 ，te=0 的 等 效 时 间 线 称 为 参考 时 间 线 。 在 瞬时 压缩 
线 上 的 应 力 应 变 关 系 为 

Ce 

E2 一 2 + 对 | (5) 


ON 


二 


u 


在 参考 时 间 线 上 的 应 力 应 变 关系 为 
sp = EP :| 到 | (6) 


oO; 0 


在 某 一 恒定 的 有 效应 力作 用 下 ， 蠕 变 可 以 表达 为 
E2 = 多 a “| (7) 
大 


to 
式 中 yy /V 为 表征 土 的 蚂 变 能 力 的 蚂 变 参数 ， 不 同 于 
里 变 系数 Ce。 里 变 率 为 


(8) 


当 t= 时 ，22? 趋 近 于 无 穷 小 ， 不 产生 蚂 变 ， 称 为 极 
限时 间 线 。 如 图 2 中 点 i 与 点 i” 所 示 , 超 固 结 点 等 效 
时 间 较 大 ， 因 而 蠕 变 率 较 小 ， 正 常 固 结 点 等 效 时 间 较 
小 ， 因 而 蜂 变 率 较 大 。 任 意 一 点 的 应 变 可 表示 为 

| [a ea (9) 


Oz0 V 10 


旧式 〈9) 可 得 


本 


1 
nom- | (10) 
wv lo 


对 于 任 一 应 力 应 变 状态 点 ， 有 


xdo+exnd (1D 
Vo. 


de, =ds +de®? = 


将 式 (8)、(10) 代入 式 (11)， 得 到 


1 
。 1 V i 
Ee Se + exp 人 a < (12) 
Vo VW Ww oo 


式 〈12) 即 是 描述 士 的 一 维 弹 黏 塑性 本 构 模 型 的 微分 
方程 ,应 变 率 分 为 与 时 间 无 关 的 弹性 率 和 与 时 间 有 关 
的 黏 性 率 。 理 论 上 ， 其 所 反映 的 应 力 一 应 力 率 一 应 变 
一 应 变 率 关系 适用 于 任何 加 载 形 式 。 殷 建华 采用 该 模 
型 分 别 给 出 了 等 应 变 率 试 验 、 等 应 力 率 试验 、 应 力 松 
弛 试验 等 试验 条 件 下 应 变 的 计算 方法 051。 


“、(o' ,£2) Ino. 
. 20?“:0 3 


/ | 
(G1,82 | 
| 


瞬时 压缩 时 间 线 
(CE 
极限 时 间 线 一 w\ 


等 效 时 间 线 一 < 


2 应 力 应 变 对 数 关系 [1 芭 
Fig. 2. Relationship of vertical strain versus In(vertical effective 


stress). 


1.3 考虑 流 变性 的 非 线 性 固 结 模型 
使 用 有 效应 力 原理 af +x = 可 以 建立 式 (4) 的 
超 静 孔隙 水 压力 和 式 〈12) 中 有 效应 力 之 间 的 关系 ， 
使 一 维 弹 妖 塑 性 本 构 关系 与 土 中 孔隙 水 压力 的 非 线 性 
消散 过 程 焰 合 。 进 而 得 到 同时 包含 se.(z,z)、ulz?) 的 偏 
微分 方程 组 , 即 为 土 的 一 维 弹 黏 塑性 非 线性 固 结 模型 。 
对 方程 组 中 的 常量 赋值 并 给 出 初始 条 件 、 边 界 条 件 ， 
即 可 得 出 土 体 在 固 结 过 程 中 应 变 、 超 静 孔 阶 水 压力 与 
深度 z、 时 间 1 的 函数 关系 。 
一 维 EVP 非 线性 固 结 模型 
条 件 为 


微分 方程 及 初始 、 边 界 


1 

-人 | 人 ea (13) 
u(0,1)=0 
u(H,t)=0 
u(z,0)= 超 静 孔 压 初始 条 件 
s.(z,0)= 应变 初始 条 件 
式 中 : 瑟 为 承载 土 层 的 厚度 。 

一 维 非 线性 固 结 与 一 维 EVP 本 构 模 型 耦合 分 析 
得 出 的 应 变 是 土 体 深度 和 时 间 的 函数 ， 代 表 某 一 时 刻 
土 体 中 某 一 点 处 的 应 变 值 ， 而 在 实际 工程 与 土工 试验 
中 ， 往 往 更 加 关注 土 体 的 总 应 变 。 本 文采 用 式 (14) 
计算 承载 土 体 的 总 应 变 。 


a = |, i) (14) 


2 一 维 循 环 荷 载 国 结 试验 


2.1 试验 简介 
杭州 地 区 底层 履 盖 大 量 的 深厚 软 昔 土 ， 长 期 地 铁 
等 循环 荷载 作用 下 引起 的 沉降 问题 不 容 忽视 。 本 文 为 


验证 一 维 弹 黏 塑性 非 线性 固 结 模型 ， 针 对 杭州 地 区 软 
黏 士 进行 了 一 系列 的 一 维 循环 荷载 固 结 试验 。 


试验 装置 采用 的 动力 循环 蠕 变 试验 仪 。 其 主要 衣 
件 包含 固 结 容器 、 加 载 系统 、 控 制 系统 、 量 测 系统 和 
数据 采集 系统 。 采 用 标准 的 固 结 容 器 ， 其 尺寸 为 高 
20mm， 直 径 61.8mm。 加 载 系统 通过 无 电 刷 直流 伺服 
马达 控制 器 和 高 精度 位 移 传 感 器 形成 闭环 控制 。 整 个 
系统 由 电脑 进行 自动 控制 ， 可 施加 余弦 荷载 、 三 角形 
筒 载 、 梯 形 荷 载 等 。 本 文 试验 用 土 为 杭州 软 条 土 的 饱 
和 重 塑 土 ， 含 水 量 高 ， 流 变 特性 显著 ， 物 理 参数 见 表 
1 。 


在 可 选 荷载 形式 中 ， 余 弦 荷 载 相对 于 其 他 形式 的 
循环 荷载 ， 具 有 高 阶 可 导 性 ， 便 于 模型 计算 。 本 试验 
中 土 样 预 固 结 6h 后 ,连续 施加 20 个 周期 的 余弦 荷载 ， 
其 中 试验 1 在 土 样 预 固 结 后 连续 作用 4 种 不 同形 式 的 
荷载 ， 有 具体 荷载 形式 见 表 2。 


表 1 杭州 软 黏 士 的 物理 参数 
Table 1 Physical parameters of soft clay in Hangzhou. 
液 限 塑 限 塑性 指数 比重 
1 WP Lh Gs 


40.3% 24.2% 16.1% 2.73 


1.32 Mg/ms 


表 2 荷载 形式 与 土 样 初始 状态 


Table 2 Load form and soil sample initial state. 


预 固 结 压 
人 答 载 形式 初始 孔 院 
试验 组 号 力 Pe 
Peyc(kPa) 比 e0 
(kPa) 
100+10cos(0.12xt) 
100+20cos(0.12xt) 

1 100 1.097 
100+30cos(0.12xt) 
100+40cos(0.12xt) 

2 100 100+10cos(0.24rrb) 1.187 

3 100 100+10cos(0.48rb) 1.325 

4 200 200+20cos(0.12xt) 1.040 

5 300 300+30cos(0.12xt) 1.126 


2.2 ”模型 参数 的 确定 

应 当 认 识 到 ， 土 的 固 结 模型 的 提出 是 为 了 解决 实 
际 工程 的 需要 ， 模 型 的 准确 性 和 使 用 方便 性 都 应 得 到 
足够 的 重视 。 为 了 使 本 文 提出 的 一 维 弹 攻 塑性 非 线性 
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3 e-lgh 拟 合 曲 线 
Fig. 3. Relationship of void ratio versus lg(Permeability 


coefficient). 


o%、e% 是 参考 时 间 线 上 的 参考 点 , 以 确定 出 入 线 
的 位 置 ， 本 文中 土 体 的 基本 力学 参数 由 标准 固 结 试验 
测 得 , 因而 可 由 单 级 加 载 1day 的 一 维 固 结 试验 进行 测 


回 结 模型 能 与 土工 试验 较 好 地 结合 ， 便 于 应 用 。 本 文 


部 分 参考 括 建 华 提 出 的 假设 B 简化 算法 P3, 推荐 采用 
以 下 确定 参数 的 方法 。 
Kk、4 可 通过 标准 固 结 试验 确定 ， 满 足以 下 关系 : 
K C 
一 三 一 一 15 
V VIin1l0 3 
1 C 
一 二 一 一 人 16 
V VIinl0 oe) 


其 中 , ce 是 回 弹指 数 。kw 是 有 效应 力 等 于 0o0 时 ， 土 体 
的 竖 向 渗透 系数 。o 取 每 组 试验 的 动态 偏 移 量 , hw 可 
由 GDS 先进 固 结 试验 系统 测量 R89, 若 试验 条 件 受 限 ， 

可 采用 e-lghk, 拟 合 直 线 求 得 不 同 孔隙 比 下 的 渗透 系数 。 
曾 巧 玲 通 过 对 7 种 固 结 系数 确定 方法 的 计算 结果 与 试 
验 结果 进行 对 比 ， 指 出 三 点 法 具有 精度 较 高 ， 计 算 简 
便 等 优点 所 ,文献 [26] 指 出 使 用 三 点 法 计算 渗透 系数 
误差 较 小 。 因 此 ， 本 文 对 三 组 标准 固 结 试验 的 试验 结 
果 采 用 三 点 法 计算 固 结 系 数 ， 并 由 式 h=Cyx mvxyw 计 
算出 各 孔隙 比 下 的 渗透 系数 。 其 中 体积 压缩 系数 wz 由 
各 级 固 结 试验 的 初始 与 结束 时 刻 应 变 的 差 值 得 到 ， 在 
同一 组 固 结 试验 中 不 取 和 定 值 。 拟 合 出 的 e-lgh 直线 见 
图 3， 拟 合 直 线 的 斜率 即 为 渗透 指数 ce， 因而 有 
cF=0.2694 。 本 文 试 验 中 压缩 指数 cc=0.2466， 则 
cocr=0.915。 已 知 土 在 特定 应 力 应 变 状态 下 的 孔隙 比 ， 
可 在 e-lgh 拟 合 直 线 上 插值 ， 得 到 ko。 


量 :将 制备 好 的 土 样 放 入 固 结 容器 ,施加 某 一 静 载 o'， 
且 o', 大 于 先期 回 结 压力 , 测 出 1440 min 时 的 应 变 值 ， 
即 为 sw 。 

ts 取 固 结 试验 单 级 加 载 时 间 , 本 文中 t=1440 min。 

% 是 表征 土 的 流 变 能 力 的 蠕 变 参数 ， 在 一 维 EVP 
本 构 模型 中 ， 其 值 无 法 直接 测量 ， 可 通过 与 固 结 试验 
结果 拟 合 得 到 , 详 见 文献 [19]。 在 准确 性 要 求 不 高 时 ， 
可 采用 式 (17)〉 近似 计算 。 

Coe 
i 840 

本 文 试 验 中 ， 获 得 的 杭州 软 黏 士 的 参数 x 为 
0.0153， 4 为 0.1071， 里 变 参 数 w 为 0.00174。 模 型 计 
算 中 用 的 其 它 各 项 参数 见 表 3。 


2.3 ”初始 条 件 的 确定 

由 土 样 的 标准 固 结 试验 结果 可 知 ， 预 固 结 6h 后 ， 
主 固 结 已 完成 ， 因 而 第 一 级 循环 加 载 的 初始 应 变 采用 
常数 ， 数 值 上 等 于 预 固 结 结束 时 的 应 变 值 。 由 于 试验 
装置 的 轴 力 控制 系统 只 能 施加 连续 荷载 ， 在 预 固 结 结 
束 后 ， 应 力 会 从 预 固 结 压力 以 较 快 的 速度 连续 增加 到 
循环 荷载 的 初 值 ， 从 而 造成 相位 偏差 。 为 了 使 相位 偏 
差 的 影响 最 小 化 ， 计 算 时 采用 原 荷载 形式 左 移 1/8 个 
周期 后 的 荷载 形式 ， 见 图 4。 

对 于 试验 1， 第 二 阶段 及 之 后 的 各 个 阶段 的 初始 
应 变 与 初始 孔 压 可 采用 前 一 阶段 结束 时 刻 的 计算 结果 。 
使 用 Mathematica 求解 式 〈13) 可 得 出 第 一 阶段 的 应 
力 一 深度 一 时 间 和 应 变 一 深度 一 时 间 数 值 结果 。 导 出 
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第 一 阶段 结束 时 刻 的 应 力 一 深度 和 应 变 一 深度 的 数据 ， 
并 拟 合 出 多 项 式 趋势 线 ， 作 为 下 一 阶段 的 初始 条 件 。 


第 二 至 第 四 阶段 的 初始 超 静 孔 际 水 压力 拟 合 结果 见 图 


3 。 


由 于 本 文中 的 模型 适用 于 任何 单 向 加 载 形 式 ， 此 
初始 条 件 确定 方法 在 其 他 形式 的 多 级 加 载 固 结 试验 中 


同样 适用 。 


表 3 计算 模型 中 各 项 参数 值 


Table 3 The parameter values in the calculation model. 


3 ”模型 验证 


为 方便 模型 的 计算 ， 绘 制 了 使 用 
件数 值 求解 式 〈13) 的 程序 术 


开始 
定义 各 模型 参数 变 
量 ， 并 赋值 
于 
编写 偏 微分 方程 组 ， 给 出 
初始 、 边 界 条 件 


学 


数值 求解 深度 、 时 间 相 关 


的 偏 微分 方程 组 


| 
Y 


积分 得 到 沿 土 层 厚度 的 
应 变 ， 获 得 总 应 变 - 时 间 


臣 图 ， 见 图 6。 


更 新 荷载 参数 、 渗 透 
系数 和 初始 条 件 


人 A 
更 新 下 级 加 载 的 初 
始 条 件 


术 


Mathematica 软 


试验 组 kvo az0 
号 (m/min) (kPa) 
0.62e-8 
1 100 0.095 
0.57e-8 
0.55e-8 
2 1.31e-8 100 0.092 
3 2.80e-8 100 0.114 
4 0.36e-8 100 0.065 
5 0.31e-8 100 0.081 
| 一 计算 各 
一 “试验 荷载 


计算 荷载 与 试验 
荷载 的 差异 


| , 1 1 
360 362 364 366 368 
时 间 (min) 


图 4 相位 差 调整 图 示 
Fig. 4. Adjustment ofphase difference. 


EE. 
. 第 二 级 初始 孔 压 ; 

gj Sa =0.121x°-2.193x+0.048 (RD 
© -12 一 纪 Te 
名 第 三 级 初始 孔 压 : 
尺 SY 3=0.00017x-0.00613€+0.031x- 2 
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370 


5 各 阶段 初始 孔 压 拟 合 结果 


Fig. $. Fitting results of initial excess pore water pressure in each 


Stage. 


导出 应 变 、 时 间 关 系 的 导出 结束 时 刻 的 应 变 - 深 
数据 度 和 和 孔 压 -深度 数据 
y 和 A 
继续 加 载 ? 是 
否 
4 
结束 


6 模型 求解 流程 图 


Fig. 6. Flow chart of model solving. 


3.1 和 孔 压 分 析 


图 7 为 试验 1 第 一 级 循环 荷载 作用 下 前 7 个 周期 
的 超 静 孔隙 水 压力 模型 解 。 从 图 中 结果 可 以 看 出 ， 在 


幅 越 大 ， 这 与 土 的 固 结 理论 相 一 致 。 


土 样 中心 点 处 的 超 静 孔 压 波动 


循环 荷载 作用 下 ， 了 筷 际 水 压力 随 着 时 间 的 增长 在 正 负 
孔 压 之 间 振 动 并 趋 于 稳定 。 在 土 体 的 排水 边界 ， 孔 阶 
水 压力 为 0， 在 土 体 1/2 深度 处 《〈 排 水 距离 最 远 )， 孔 
压 振 幅 最 大 。 在 土 体 的 竖 向 上 ， 距 离 排 水 边界 越 远 的 
部 位 孔隙 水 的 排出 与 进入 越 缓 慢 ， 因 而 超 静 孔 压 的 振 


昌 线 如 图 8 所 示 。 


随 着 荷载 振幅 的 增加 ， 超 静 孔 隙 水 压力 的 振幅 也 相应 
增加 ， 谷 值 绝对 值 与 荷载 振幅 基本 相等 ， 而 峰值 则 逐 
渐 超 过 荷载 振幅 。 由 于 黏 塑性 的 存在 ， 在 循环 荷载 作 
用 下 ， 土 体 不 但 产生 弹性 应 变 波动 ， 塑 性 应 变 也 会 随 


着 作用 时 间 而 不 断 增 大 。 同 时 ， 


在 土 样 中 心 处 由 于 和 孔 


隙 水 不 易 排 出 ， 有 效应 力 波动 很 小 ， 土 颗粒 骨架 会 不 
断 地 产生 里 变 ， 使 有 效应 力 发 生 微小 的 松弛 现象 ， 又 


由 于 竖 向 力 平衡 ， 而 使 得 孔 际 水 压力 的 最 大 值 超过 蓓 


载 最 大 振幅 。 此 现象 可 以 从 孔 隐 水 压力 与 循环 荷载 的 


相位 关系 中 “如 图 9) 进一步 的 验证 。 从 图 9 可 以 看 


出 ， 土 样 中心 处 超 静 孔隙 水 压力 与 循环 荷载 存在 一 定 
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的 相位 偏差 。 蓓 载 增 大 会 引起 孔隙 水 压力 的 增 大 ， 同 


形 特 征 ， 本 文采 用 在 不 同 频 率 、 不 同 振幅 、 不 同 循环 


时 ， 孔 隙 水 的 排除 会 造成 孔隙 水 压力 的 消散 。 由 于 和 荷 
载 的 变化 与 孔隙 水 压力 的 消散 均 是 动态 过 程 ， 最 大 孔 
隙 水 压力 并 非 在 最 大 荷载 作用 时 刻 产 生 ， 而 是 出 现在 


荷载 增 大 速率 与 孔 际 水 压 消散 速率 的 差 值 最 大 的 时 刻 。 


同 理 ， 最 大 负 孔 压 出 现在 荷载 减 小 速率 与 负 孔 压 消散 
速率 的 差 值 最 大 的 时 刻 。 
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图 7 试验 1 第 一 级 超 静 孔隙 水 压力 一 深度 一 时 间 三 维 图 像 
Fig. 7 Three-dimensional image of excess pore water pressure- 


depth-time of the first stage in experiment 1. 


超 静 孔隙 水 压力 (kPa) 


8 试验 1 土 样 中 心 点 处 超 静 孔隙 水 压力 一 时 间 图 像 
Fig. 8. Relationship of excess pore water pressure versus time at 


the center point of soil sample in experiment 1. 


i A 


应 力 (kPa) 


”一 起 前 孔 阶 水 压力 
; - - - 循环 荷载 


on 
时 间 (min) 
9 超 静 孔隙 水 压力 与 循环 荷载 相位 差 
Fig. 9. The phase difference between the excess pore water 


pressure and the cyclic load. 


3.2 ”模型 验证 
合理 的 模型 应 能 反映 材料 在 不 同 受 力 情 况 下 的 变 


荷载 偏 移 量 作 用 下 得 到 的 试验 结果 与 模型 计算 结果 进 
行 对 比 。 图 10 为 试验 1 第 一 级 荷载 作用 下 前 7 个 周 
期 的 应 变 一 深度 一 时 间 三 维 图 像 。 在 循环 荷载 作用 下 ， 
应 变 值 波动 式 增 长 ， 靠 近 排 水 边界 处 ， 孔 隙 水 压力 的 
振幅 较 小 ， 则 有 效应 力 的 均值 和 振幅 均 较 大 ， 因 而 应 
变 波 动 剧 烈 ， 均 值 较 大 。 在 土 样 中 心 处 ， 超 静 孔 隙 水 
压力 的 振幅 与 荷载 接近 ， 存 在 较 小 的 相位 差 ， 有 效应 
力 振幅 较 小 ， 因 而 应 变 振 幅 较 小 。 


10 应 变 一 深度 一 时 间 三 维 图 像 


Fig. 10. Three-dimensional image of strain-depth-time. 


图 11~13 分 别 给 出 了 不 同形 式 的 余弦 荷载 作用 下 
由 试验 和 模型 求解 得 出 的 总 应 变 一 时 间 曲 线 以 及 两 者 
之 间 的 最 大 偏差 。 从 试验 结果 可 以 看 出 ， 随 着 加 载 时 
间 的 增 大 ， 总 应 变 值 波动 增长 。 在 同一 种 形式 的 荷载 
作用 下 ,弹性 应 变 振幅 保持 不 变 , 塑性 应 变 持续 增长 ， 
但 增长 率 逐 渐 减 小 ， 出 现 应 变 硬化 。 

如 图 4 所 示 ， 在 试验 仪器 循环 加 载 的 初始 时 刻 ， 
计算 所 采用 的 荷载 形式 与 试验 中 加 载 的 形式 有 一 定 的 
差别 ， 前 者 是 不 连续 的 ， 而 后 者 是 连续 的 。 因 此 ， 在 
第 一 个 波段 的 计算 结果 与 试验 结果 均 有 较 大 差别 。 在 
多 级 加 载 试验 中 ， 由 当前 循环 荷载 转换 为 下 一 级 循环 
荷载 时 仍 会 由 于 试验 仪器 只 能 连续 变 荷 而 带 来 误差 。 
每 一 级 循环 荷载 施加 结束 后 ， 压 力 会 迅速 减 小 为 动态 
偏 移 量 ， 再 以 较 快 的 速率 增 大 到 下 一 级 荷载 初 值 ， 因 
而 试验 结果 的 最 后 一 个 波形 明显 不 同 于 其 他 波形 。 
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11.15 试验 结果 
划 计算 结果 
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11 不 同 频率 的 循环 荷载 固 结 试验 结果 : (a) 广 0.001 Hz; (b)/ 
=0.002 Hz; and (c)/ = 0.004 Hz 


Fig. 11. Test results of cyclic loading at different frequencies: (a) 


大 0.001 Hz; (b)f=0.002 Hz; and (ce)f = 0.004 Hz 
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12 不 同 振幅 的 循环 荷载 固 结 试验 结果 : (a) 4= 士 10 kPa; (b) 


A=+20 kPa; (c)A= +30 kPa; and (d)A= +40kPa 
Fig. 12. Test results of cyclic loading of different amplitudes: (a) A= 
士 10 kPa; (b) A=+20 kPa; (c)4= +30 kPa; and (d)A= +40kPa 
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13 不 同 动态 偏 移 量 的 循环 荷载 国 结 试验 结果 : (a) Pe= 100 
kPa; (b) Pa =200 kPa; and (c) Pa = 300kPa 
Fig. 13. Test results of cyclic loading of different dynamic offset: 
(a) Pa= 100 kPa; (b) Pa =200 kPa; and (c) Pa = 300kPa 


如 图 12 所 示 ， 试 验 1 在 同一 土 样 上 共 施 加 了 四 
级 循环 荷载 。 各 级 荷载 频率 均 为 0.001Hz, 振幅 依次 增 
大 , 分 别 为 : 士 10kPa、 士 20kPa、 士 30kPa、 士 40kPa。 
从 试验 结果 可 以 看 出 ， 由 第 一 级 至 第 四 级 ， 土 的 塑性 
应 变 增长 率 逐 渐 减 小 。 尽 管 荷载 振幅 依次 增 大 ， 但 随 
着 循环 荷载 作用 时 间 的 增加 和 士 体 的 密实 ， 依 旧 产 生 
了 较为 显著 的 应 变 硬 化 现象 。 计 算 结果 则 表现 出 更 为 
显著 的 应 变 硬化 ,计算 出 的 应 变 值 由 比试 验 结果 偏 大 ， 
在 第 四 级 加 载 阶段 变 为 偏 小 于 试验 结果 ， 且 塑性 应 变 
增长 速率 明显 小 于 试验 结果 。 本 文 研究 发 现 ， 模 型 中 
总 应 变 增长 的 形态 主要 受到 蠕 变 参数 y 值 的 影响 。 由 
试验 结果 与 计算 结果 的 对 比 可 以 大 致 判断 出 循环 荷载 
作用 对 土 样 的 蠕 变 特 性 造成 了 一 定 的 影响 。 此 外 ， 一 
些 学 者 通过 试验 研究 得 出 了 循环 荷载 作用 会 造成 士 的 
模 量 弱化 的 结论 [27]， 文 献 [1H] 中 指出 ，w 值 的 大 小 会 
受到 荷载 作用 时 间 、 应 变 大 小 的 影响 。 

图 11、12、13 中 的 对 比 结果 显示 ， 在 不 同 频率 、 
不 同 振幅 、 不 同 动态 偏 移 量 的 试验 条 件 下 ， 模 型 计算 


结果 与 试验 结果 吻合 较 好 ， 能 够 较 好 的 反映 各 种 荷载 
形式 作用 下 土 的 力学 行为 ， 验 证 了 模型 的 合理 性 。 


3.3 ”模型 的 讨论 
若 不 考虑 土 的 流 变 性 , 可 运用 传统 的 e - 1go' 曲 


线 求解 一 维 非 线性 固 结 的 变形 问题 在 塑性 压缩 阶段 ， 
压缩 系数 m, =c。 /+e,)In100']， 可 得 2 


O05 本 Ce Ou dp 18 

Ot (I+e,)Inl0o’\ or © 
在 回 弹 阶段 ，m, =c,/[4+e,)n10o']， 可 得 

加 0c Ce Ou dp 19 

ot na 人 


将 式 (18)、(19) 代入 式 (4) 的 左 侧 可 以 分 别 解 出 循 


流 变 特性 ， 无 论 荷 载 大 小 ， 都 会 发 生 么 塑性 变形 ， 只 
是 肾 变 率 不 同 。 随 着 时 间 的 增 大 ， 本 文 模型 计算 出 的 
总 应 变 持续 增 大 ， 且 与 试验 结果 相 接近 ， 而 不 考虑 流 
变性 的 计算 结果 随 着 时 间 的 增长 越发 小 于 实际 应 变 值 。 
在 长 期 循环 荷载 作用 下 ， 软 黏土 的 流 变 会 更 加 显 蓝 ， 

成 为 预测 地 基 沉 降 时 不 可 缺少 的 一 个 考虑 因素 。 
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4 结论 

本 文 将 一 维 非 线性 固 结 理论 与 一 维 弹 秋 塑性 本 
构 模型 耦合 分 析 ， 得 出 了 具有 普遍 适用 性 的 一 维 弹 秋 
塑性 非 线性 固 结 模型 。 通 过 进行 一 系列 不 同 频率 、 不 
同 振幅 、 不 同 动态 偏 移 量 的 循环 荷载 固 结 试验 ， 对 比 
分 析 了 试验 结果 和 模型 计算 结果 。 对 比 结果 表明 ， 在 


环 荷载 作用 下 土 体 在 压缩 和 回 弹 阶 段 的 超 静 孔隙 水 压 
力 。 根 据 有 效应 力 原理 ， 可 以 得 到 土 体内 部 有 效应 力 
的 变化 。 基 于 土 体 的 压缩 、 回 弹 曲 线 ， 如 下 : 


| (压缩 曲线 ) 
co Go (20) 


2 ( 回 弹 曲线 ) 
lt+e, Na 


€=6,+ 
1+ 


可 以 得 到 应 变 在 不 同 深度 随时 间 的 变化 特征 ， 进 而 由 
式 (14) 求 得 在 不 同时 刻 土 体 的 总 应 变 。 用 上 述 忽 略 
流 变 的 传统 方法 计算 本 文 试验 4 的 应 变 发 展 ， 计 算 结 
果 见 图 14。 


x 


10.8 -5 


2 ”试验 结果 
一 一 一 本文 模型 计算 结果 
-不 考虑 流 变 的 计算 结果 


应 变 (%) 


nl 
0 50 100 150 20 250 300 350 
时 间 (min) 
14 本 文 模型 与 传统 模型 计算 结果 对 比 
Fig. 14. The comparison between the calculation results of the 


models in this paper and ignoring rheology is made. 


在 循环 荷载 作用 下 ， 由 不 考虑 流 变 性 的 一 维 非 线 
性 固 结 理论 计算 结果 显示 : 土 体 在 加 载 初始 阶段 因 有 
效应 力 的 增 大 而 产生 较 大 的 塑性 变形 ， 之 后 则 以 恒定 
的 振幅 持续 波动 ， 且 均值 不 变 ， 仅 能 反映 弹性 变形 的 
振动 特点 ， 而 不 能 反映 黏 塑性 变形 的 加 剧 。 而 本 文 提 
出 的 一 维 非 线性 弹 秋 塑性 固 结 模型 由 于 考虑 了 土 体 的 


不 同 的 荷载 形式 下 ， 模 型 的 计算 结果 均 可 以 较为 合理 
的 预测 土 体 在 循环 荷载 作用 下 的 应 变 发 展 。 主 要 的 结 
论 如 下 : 
(1) 耦合 得 到 的 一 维 弹 秋 塑 性 非 线 性 固 结 模型 
能 够 表征 土 体 在 固 结 过 程 中 的 超 静 孔隙 水 压力 一 深 
度 一 时 间 关 系 和 应 变 一 深度 一 时 间 关 系 。 确 定 各 参数 
值 、 初 始 条 件 和 边界 条 件 后 可 以 用 此 模型 对 土 体 任 一 
时 刻 任 一 深度 处 的 超 静 孔 压 和 应 变 求解 。 
(2) 为 了 使 模型 更 好 的 与 土工 试验 相 结合 , 给 出 
了 即便 于 获取 又 能 较 准 确 地 反映 土 体 变形 特性 的 模 
型 参数 确定 方法 ， 并 且 给 出 了 多 级 单 向 加 载 试验 的 初 
始 条 件 确 定 方法 。 
(3) 在 循环 荷载 作用 下 ， 超 静 孔 隙 水 压力 、 有 
效应 力 、 应 变 均 表 现 为 波动 式 增长 。 在 土 体 中 排 
水 距离 较 远 的 部 位 ， 循 环 荷载 与 超 静 孔隙 水 压力 会 因 
为 固 结 作用 而 产生 一 定 的 相位 差 ， 且 有 效应 力 的 振幅 
和 应 变 振幅 随 排水 距离 的 增 大 而 减 小 。 随 着 循环 荷载 


作用 时 间 的 增 大 ， 土 体会 产生 应 变 硬化 现象 ， 总 应 变 
增长 趋势 变 缓 。 
(4) 通过 对 比 五 组 一 维 循环 荷载 固 结 试验 的 试 


验 结 果 与 模型 的 计算 结果 ， 表 明 由 本 文中 模型 计算 出 
的 总 应 变 可 以 较为 合理 的 预测 土 体 在 循环 荷载 作用 
下 的 变形 过 程 。 
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